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INTRODUCTION 

During  the  past  year  two  of  our  students  completed  research  projects. 

Edward  Josberger  finished  his  Ph.D.  thesis  on  the  ablation  of  vertical  Ice 
walls  submerged  In  warm  sea  water,  and  Donald  Perovich  finished  his  M.S. 
thesis  on  the  optical  properties  of  young  sea  Ice.  Scientific  reports  on 
both  topics  have  been  prepared  and  distributed.  Calculations  of  the  large- 
scale  heat  and  mass  balance  over  the  AIDJEX  array  have  been  carried  out  and 
compared  with  values  in  the  central  part  of  the  Basin  as  Inferred  from 
earlier  drifting  station  data.  The  sensitivity  of  regional  fluxes  to 
processes  associated  with  Ice  deformation  was  also  Investigated.  Laboratory 
measurements  of  light  scattering  by  sea  Ice  have  been  concluded.  Results 
Indicate  that  the  scattering  coefficient  depends  on  both  temperature  and 
growth  rate,  while  the  phase  function  was  found  to  be  unexpectedly  insensitive 
to  these  quantities.  Analysis  of  laboratory  and  field  work  on  the  Interaction 
of  waves  and  grease  Ice  has  been  completed,  leading  to  the  formulation  of  an 
empirical  equation  for  wave  damping  by  grease  Ice.  A  theoretical  model  of 
solar  heating  In  leads  Is  being  developed.  Preliminary  work  suggests  that  the 
rate  of  lateral  melting  responds  directly  to  changes  In  lead  width.  If  the  lead 
Is  narrow;  however.  Increases  In  width  above  about  100  m  have  little  effect 
because  the  greater  heat  Input  Is  balanced  by  greater  heat  loss  to  the  atmosphere. 

A  spring  cruise  on  the  NOAA  ship  SURVEYOR  was  carried  out  In  the  Bering 
Sea.  Wave  attenuation  at  the  Ice  edge  and  the  relative  motions  of  Ice  bands 
which  break  free  from  the  pack  were  Investigated.  The  new  Infrared  scanning 
photometer  was  completed  In  April  and  successfully  field  tested  at  Barrow 
during  May  and  June.  Changes  In  spectral  albedos  and  Incident  Irradlance 
as  a  function  of  cloudiness  were  recorded  through  the  early  melt  stages. 
Experimental  methods  for  measuring  lateral  ablation  In  leads  were  also  tested. 
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Scientific  results  from  previous  work  Include  papers  on  the  effects  of 
Ice  thickness  on  the  absorption  of  solar  radiation,  brine  drainage  and 
convection  In  young  Ice,  the  melting  of  vertical  Ice  walls  In  warm  salty 
water,  the  optical  properties  of  young  Ice,  Ice  edge  phenomena,  and  wave 
damping  by  grease  Ice. 


DYNAMIC  AND  THERMODYNAMIC  MODELING 

Heat  and  Mass  Balance  Studies.  Our  goals  here  are:  (i)  to  understand 
large-scale  Ice  production  and  surface  heat  exchange  over  the  Arctic  Ocean, 
and  (11)  to  determine  the  sensitivity  of  these  quantities  to  processes 
associated  with  the  deformation  of  the  Ice  pack.  During  the  past  year  we 
have  been  engaged  In  a  comprehensive  study  of  this  problem,  utilizing  all 
available  data.  The  basic  tool  for  this  work  Is  our  regional  averaging 
mode!  which  requires:  (1)  Ice  strain  data  to  calculate  the  thickness 
distribution  and  (11)  Incident  energy  fluxes  to  calculate  Ice  growth  and  the 
surface  heat  balance  for  each  thickness  of  Ice.  Since  the  early  calculations, 
the  model  has  been  Improved  through  the  addition  of  Improved  numerical 
techniques  and  more  realistic  ridging  parameters.  Testing  these  Improvements 
with  the  T-3/ARLIS  II/  NP-10  strain  data  produced  a  substantial  Increase  In 
the  predicted  fluxes.  Regional  fluxes  have  also  been  calculated  using  high 
quality  AIDJEX  strain  data  and  were  found  to  be  significantly  different  than 
those  from  the  T-3  triangle.  For  example,  annual  totals  of  sensible  heat 
loss  were  two  and  a  half  times  larger  over  the  AIDJEX  manned  array,  while  Ice 
production  was  35%  greater.  The  latent  heat  flux,  on  the  other  hand, 
differed  by  only  8%.  AIDJEX  strains  tended  to  put  more  Ice  In  the  20-40  cm 
thickness  category  and  less  Into  the  40-80  cm  category,  causing  area-weighted 
turbulent  flux  values  to  peak  In  the  20-40  cm  range  rather  than  the  40-80  cm 
range  found  with  the  T-3  data.  It  Is  not  clear  whether  these  large  differences 
result  from  normal  year-to-year  changes  In  Ice  movement,  differences  In  Ice 
flow  between  the  Beaufort  Sea  and  the  Central  Arctic,  Improvements  In 
measurement  accuracy,  differences  In  the  scale  of  the  two  sets  of  strain 
measurements,  or  some  combination  of  the  above.  Upcoming  data  from  the  FGGE 
buoy  array  should  be  helpful  In  resolving  this  question. 


In  the  second  phase  of  this  work,  we  have  carried  out  over  a  dozen 
additional  simulations  in  order  to  learn  more  about  how  the  above  results 
are  affected  by  strain  field  characteristics  and  by  assumptions  regarding 
the  ridging  parameters.  We  found  that  setting  the  shear  to  zero  and  leaving 
the  divergence  unchanged  in  either  the  T>3  or  AIDJEX  data  caused  about  a 
50%  reduction  in  the  amount  of  open  water  during  the  winter  months.  Thus, 
although  shear  dominates  the  strain  field  of  the  ice  pack,  such  motions  are 
responsible  for  no  more  than  about  half  of  the  winter  open  water  production. 
During  the  sunmer,  divergence  and  lateral  melting  completely  determine  the 
amount  of  open  water.  Despite  decreases  in  the  amount  of  open  water  and 
young  ice  when  shear  was  neglected,  there  was  only  about  a  1D%  decrease  in 
the  winter  rates  of  ice  production.  This  was  because  of  substantially  larger 
amounts  of  thick  first-year  ice  which  partially  compensated  for  the  ice 
production  lost  in  the  thinner  categories.  The  sensible  heat  flux,  on  the 
other  hand,  responded  strongly  to  the  decrease  in  thin  ice;  the  magnitude  of 
the  change  depended  on  season,  but  for  the  year  there  was  over  a  6D%  decrease 
in  the  sensible  heat  input  to  the  atmospheric  boundary  layer.  Dther 
calculations  indicated  that  the  regional  heat  balance  is  quite  sensitive  to 
assumptions  regarding  which  thicknesses  of  ice  participate  in  ridging,  but  is 
insensitive  to  the  way  in  which  this  ice  is  redistributed.  An  unexpected 
result  was  that  the  large-scale  heat  fluxes  generated  from  daily  AIDJEX  strain 
data  were  essentially  unaffected  when  recalculated  with  strains  which  were 
passed  through  a  5  day  filter.  This  suggests  that  motions  with  a  period  of 
less  than  5  days  are  unimportant  in  terms  of  ice  production  and  heat  exchange. 
We  are  carrying  out  additional  calculations  to  see  if  this  is  reasonable. 


The  final  case  we  are  studying  is  a  comparison  of  the  AIDJEX  manned  array  and 
the  AIDJEX  buoy  array  to  see  if  there  are  any  consistent  differences  in 
heat  exchange  between  the  100  and  500  km  space  scales.  A  paper  describing 
the  results  of  this  work  should  be  completed  shortly. 

The  techniques  and  models  developed  for  the  above  work  have  so  far  been 
applied  only  to  the  case  of  uniform  strain  across  some  area.  It  should  also 
be  possible  to  approach  the  more  general  problem  of  a  large  region  containing 
many  strain  elements  in  much  the  same  way.  A  solution  to  this  general  problem 
is  needed  before  we  can  effectively  utilize  data  from  large  arrays  of  buoys 
We  are  therefore  working  on  the  design  of  a  program  which  will  simultaneously 
calculate  area-averaged  heat  exchange  over  many  grid  elements,  using  as  input 
only  ice  motion,  temperature,  and  cloudiness  data.  The  most  difficult  part 
of  the  problem  is  to  deduce  the  various  components  of  the  heat  and  mass  balance 
from  such  scanty  information.  Fortunately,  our  previous  work  on  improved 
ice  forecasting  equations  indicates  that  we  should  be  able  to  obtain  reasonable 
estimates  of  all  the  fluxes,  with  the  possible  exception  of  incoming  shortwave 
radiation.  We  suspect  that  parameterizations  presently  used  in  large-scale 
models  are  appropriate  for  thick  lower  latitude  clouds,  but  tend  to 
underestimate  solar  radiation  over  the  Arctic.  We  are  looking  into  this, 
possibility  and,  if  necessary,  will  develop  a  more  suitable  parameterization. 

We  are  also  trying  to  resolve  some  questions  regarding  the  advection  of  ice 
thickness  distribution  between  adjacent  grid  elements.  The  ultimate  results 
of  these  efforts  should  be  a  model  which  will  allow  us  to  determine  both 
spatial  and  temporal  changes  in  the  heat  and  mass  balance  across  the  entire 
Arctic  Basin. 


Response  of  the  Ice  to  Temperature  Changes.  Recent  concern  about  the 
increasing  levels  of  carbon  dioxide  in  the  atmosphere  and  predictions  of  a 
possible  5-10“C  increase  in  air  temperatures  in  the  polar  regions  led  us  to 
examine  the  effects  of  such  changes  on  the  ice  pack.  Our  calculations 
indicate  that  winter  temperature  anomalies  have  a  relatively  small  effect  on 
the  ice.  Summer  temperature  anomalies,  on  the  other  hand,  produced  dramatic 
changes  in  the  ice  cover.  For  example,  a  1®C  increase  in  the  10  meter  air 
temperature  during  the  melt  season  caused  roughly  a  60%  decrease  in  average 
ice  thickness  -  an  increase  of  2-3®C  caused  the  eventual  disappearance  of  the 
ice. 

The  above  results  point  up  a  weakness  implicit  in  this  approach.  While 
it  is  reasonable  to  expect  that  changes  in  climate  can  change  winter  air 
temperatures  by  several  degrees,  it  is  not  clear  that  summer  air  temperatures 
near  the  surface  of  the  melting  ice  could  be  maintained  much  above  the  freezing 
point.  Yet  it  is  the  summer  temperatures  which  seem  to  be  critical  in  deter¬ 
mining  the  long-term  response  of  the  ice.  The  question  which  now  needs  to  be 
considered  is  how  changes  in  upper  level  temperatures  (e.g.  at  the  top  of  the 
planetary  boundary  layer)  would  affect  the  ice  cover.  To  address  this 
question  we  plan  to  combine  a  suitable  boundary  layer  model  with  our  ice  growth 
model . 

Summer  Decay  of  the  Ice  Pack.  Solar  radiation  absorbed  in  the  upper 
ocean  appears  to  be  an  Important  factor  in  the  summer  retreat  of  the  ice  edge 
and  in  the  general  decay  of  the  interior  pack.  However,  precise  details  of 
the  processes  involved  are  not  known.  Lateral  melting  on  floe  edges  has 
received  some  attention  because  of  the  positive  feedback  potential  -  lateral 
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melting  causes  a  decrease  In  Ice  concentration  which  allows  more  shortwave 
radiation  absorption  In  the  water  and  a  continuing  acceleration  In  the  rate 
of  lateral  melting.  Simple  models  of  this  process  have  been  formulated  by 
Zubov  and  by  Langleben.  A  serious  drawback  with  both  these  models  Is  their 
assumption  that  all  the  absorbed  radiation  goes  Into  lateral  melting.  In 
reality  some  of  the  energy  Is  absorbed  beneath  the  Ice  where  It  contributes 
to  bottom  ablation,  while  some  fraction  of  the  remainder  Is  lost  to  the 
atmosphere  through  longwave  radiation  and  turbulent  heat  exchange. 

As  the  first  step  In  a  combined  theoretical  and  experimental  study  of 
the  problem,  we  have  developed  a  simple,  but  more  realistic,  lateral  melting 
model.  We  assume  that:  (1)  water  In  the  lead  Is  well  mixed,  (11)  energy 
absorbed  below  the  Ice  contributes  exclusively  to  melting  on  the  bottom 
rather  than  on  lead  walls,  and  (ill)  energy  deposited  In  the  layer  between 
the  surface  of  the  water  and  the  bottom  of  the  Ice  goes  Into  lateral  melting, 
warming  of  the  water,  and  heat  Input  to  the  atmospheric  boundary  layer.  An 
energy  balance  equation  which  Includes  all  the  heat  fluxes  Is  applied  at  the 
air/water  Interface,  while  lateral  melting  rates  are  determined  from  predicted 
far-fleld  temperatures  using  the  parameterization  developed  by  Ed  Josberger 
(see  "Reports  Published  and  In  Press"  section). 

The  major  difference  between  the  results  of  this  model  and  those  of 
earlier  ones  Is  In  how  the  lateral  melting  Is  affected  by  lead  width.  The 
older  models  Implicitly  assume  that  melt  rates  depend  linearly  on  lead  width; 
our  model,  however,  predicts  a  rapid  increase  in  melt  rate  only  up  to  a  lead 
width  of  about  100  meters,  after  which  melting  becomes  Increasingly  Insensitive 
to  lead  width.  Intuitively  this  seems  correct,  suggesting  that  the  Zubov/ 
Langleben  models  may  seriously  overestimate  the  speed  at  which  Ice  concentration 
decreases. 


The  problem  with  our  present  model  is  that  it  ignores  the  details  of 
heat  transport  in  the  water.  If  the  water  salinity  is  above  25  °/oo» 
radiative  heating  increases  the  stability  of  the  layer  resulting  in  a  much 
greater  eddy  diffusivity  in  the  horizontal  than  in  the  vertical.  Whether 
this  tends  to  increase  or  decrease  lateral  melting  is  not  entirely  clear. 

Another  problem  is  that,  as  lead  area  increases,  there  will  necessarily  be 
more  exchange  of  energy  between  water  in  the  lead  and  water  below  the  ice. 

Again,  we  don't  know  whether  this  exchange  will  have  a  positive  or  negative 
effect  on  the  lateral  melting.  Our  next  step  in  the  improvement  of  the  model 
will  thus  be  to  include  a  more  realistic  treatment  of  the  water  in  the 
stable  case.  This  will  be  done  through  the  introduction  of  a  two-dimensional 
heat  diffusion  equation.  Such  an  approach  has  been  tried  by  Soviet 
investigators  (e.g.  Doronin)  with  questionable  success,  probably  due  to 
inadequate  attention  to  heat  exchange  at  the  surface.  Ultimately,  we  plan  to 
treat  the  unstable  case  and  explicitly  examine  water  motion  in  leads. 

Numerous  other  mechanisms  are  potentially  important  in  the  decay  of  the 
ice  cover,  but  there  seems  to  be  no  way  to  assess  their  importance  theoretically. 
We  therefore  plan  to  carry  out  a  summer  field  program  to  determine  in  detail 
how  solar  radiation  affects  the  interaction  between  the  ice  cover  and  upper 
part  of  the  ocean.  During  the  spring  field  experiment  (see  Radiative  Transfer 
section)  several  days  were  devoted  to  testing  techniques  to  measure  lateral 
ablation  and  heat  exchange  in  summer  leads.  In  mid-June  an  ablation  array  was 
set  up  in  a  small  lead  which  had  opened  through  a  section  of  first-year  ice 
within  easy  access  from  shore.  Scuba  diving  was  carried  out  in  the  lead, 
enabling  us  to  study  the  cross-sectional  profile  of  the  ice  at  several  locations. 


The  lead  closed  before  useful  ablation  data  could  be  obtained,  however,  we 
found  that  scuba  diving  and  snorkeling  under  melt  season  conditions  was  less 
difficult  than  anticipated.  Unlike  cold  weather  diving,  a  portable  tent  was 
ample  shelter  for  the  divers,  freezeup  was  not  a  problem,  and  all  equipment 
could  easily  be  transported  on  hand  pulled  sleds.  Thus  we  feel  that  the 
only  serious  obstacle  to  a  successful  suimer  lead  experiment  Is  the  problem 
of  locating  a  suitable  site. 
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RADIATIVE  TRANSFER  STUDIES 

A  paper  entitled  "The  Effects  of  Ice  Thickness  on  the  Exchange  of  Solar 
Radiation  over  the  Polar  Oceans"  was  published  In  the  Journal  of  Glaciology. 

The  paper  describes  theoretical  calculations  of  the  absorption  and  transmission 
of  solar  radiation  by  young  growing  sea  Ice.  Two  distinct  types  of  Ice  are 
considered.  Of  primary  Importance  Is  young  growing  Ice  which  occurs  throughout 
the  Arctic  Basin  during  most  of  the  year.  Such  Ice  has  a  high  salinity.  Is 
relatively  homogeneous,  and  Is  usually  covered  with  dry  windpacked  snow.  The 
second  type,  white  Ice,  consists  of  a  highly  scattering  granu’ar  surface  layer 
above  a  thin  transition  zone  and  a  consolidated  lower  layer.  Although  white 
ice  does  not  constitute  a  major  fraction  of  the  thin  Ice  during  most  of  the 
year,  substantial  quantities  can  be  produced  by  ablation  of  first-year  Ice  at 
the  margin  of  the  Ice  pack  during  the  summer.  White  Ice  does  not  develop 
until  the  overlying  snow  layer  has  melted  away. 

Calculations  of  spectral  and  total  transmission,  absorption,  and  albedo 
were  carried  out  for  ice  from  2-80  cm  In  thickness,  with  up  to  40  cm  of  snow 
cover.  Energy  absorption  in  the  ice  was  sensitive  to  the  spectral  behavior  of 
the  albedo  over  the  entire  solar  spectrum  (400  to  3000  nm),  as  well  as  to  the 
effects  of  cloudiness  on  the  spectral  composition  of  the  Incident  irradlance. 
Since  It  was  assumed  that  clouds  remove  essentially  all  the  radiation  beyond 
1200  nm,  total  albedos  calculated  for  a  particular  Ice  type  varied  by  up  to 
30%  depending  on  the  degree  of  cloudiness.  Radiative  energy  absorption  was 
also  sensitive  to  cloud  cover  -  on  clear  days  It  was  as  much  as  50  times 
larger  In  the  upper  10  cm  of  the  Ice  than  on  cloudy  days.  Below  10  cm  most 
of  the  Infrared  wavelengths  had  been  removed  and  the  energy  absorption  was 
about  twice  as  large  on  clear  days.  In  addition,  total  transmission  to  the 
ocean  was  2  to  2.5  times  as  large  under  clear  skies. 
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Depth  profiles  of  net  irradiance  showed  that,  for  ice  thinner  than  80  cm. 
Beer's  law  provides  a  poor  description  of  the  radiation  field.  Melting  white 
ice,  for  example,  transmitted  up  to  3  times  as  much  energy  to  the  ocean  as 
Beer's  law  would  predict.  To  incorporate  the  results  of  the  present  study 
into  our  thermodynamic  ice  growth  model  without  using  an  involved  radiation 
model,  the  albedo  and  transmission  results  have  been  parameterized  as  functions 
of  snow  and  ice  thickness. 


Laboratory  Studies  of  Young  Ice.  Experiments  on  the  optical  properties 
of  thin  salty  ice  have  been  completed.  The  objective  of  the  investigation  was 
to  study  variations  in  the  albedos  and  extinction  coefficients  of  young  ice 
in  response  to  changes  in  ice  thickness,  growth  rate,  salinity,  and  temperature. 
Ten  experiments  were  carried  out.  These  included  two  runs  with  fresh  water 
at  -15°C  air  temperature,  six  runs  with  saline  water  (32  °/oo  initial  salinity) 
at  air  temperatures  ranging  from  -10°C  to  -37°C,  and  two  experiments  with  water 
of  intermediate  salinity  (16  °/oo)  at  -10®C  and  -20°C.  Each  run  lasted  from 
10  to  15  days  and  included  measurements  of  spectral  albedos  and  transmission 
profiles  as  the  ice  grew.  Concurrent  measurements  of  the  physical  state  of 
the  ice  included  vertical  temperature  and  salinity  profiles,  and  ice  thickness. 


After  the  ice  growth  was  completed  12  cm  cores  were  extracted  and  sectioned 
into  thin  vertical  and  horizontal  slabs.  Close-up  photographs  of  these  slabs 
were  taken  to  illustrate  variations  in  the  bubble  density  profile  and  to  show 
brine  drainage  features.  Additional  photographs  in  polarized  light  were  used 
to  delineate  the  crystal  structure  of  the  ice. 

In  each  of  the  fresh  ice  experiments  we  were  able  to  study  both  bubble- 
free  and  bubbly  ice.  Since  the  water  was  not  initially  saturated  with  air, 
the  ice  grew  to  about  80  mm  before  significant  bubble  nucleation  took  place. 


12 


As  the  Ice  grew  thicker,  densely  packed  bubbles  formed  in  long  vertical 
columns  rather  than  as  discrete  spheres  which  one  connonly  finds  In  salty 
Ice.  To  the  present  limits  of  accuracy,  the  optical  properties  of  bubble- 
free  fresh  Ice  are  Indistinguishable  from  those  of  pure  water.  As  the  bubbly 
layer  grew,  spectral  albedos  Increased  by  only  a  few  percent,  and  the  total 
transmission  decreased  only  slightly.  Indicating  an  extinction  coefficient  at 
500  nm  of  less  than  0.1  m~i.  Because  the  vertical  sides  of  the  columnar 
bubbles  Introduced  very  little  backscattering,  neither  the  albedo  nor  the 
attenuation  was  strongly  affected  even  though  the  vapor  volume  was  quite 
large.  If  the  same  vapor  volume  were  distributed  as  a  large  number  of  spheres 
with  about  the  same  diameter  as  the  cylinders,  backscattering  would  be  much 
greater.  Thus,  the  fresh  Ice  results  are  more  representative  of  lake  and 
river  Ice,  and  provide  only  a  lower  limit  for  the  Influence  of  bubbles  on  the 
optical  properties  of  salty  ice. 

Albedos  of  salty  Ice  Increased  much  more  rapidly  with  Ice  thickness  than 
those  of  fresh  ice.  Values  at  the  end  of  the  growth  phase  ranged  from  0.4  to 
0.9  depending  on  the  Ice  temperature  and  Initial  growth  rate.  A  significant 
temperature  dependence  was  also  found.  For  example,  the  albedo  at  500  nm  of 
25  cm  Ice  grown  at  an  air  temperature  of  -10®C  was  0.55,  but  when  the  Ice  was 
warmed  to  -2“C,  the  albedo  decreased  by  20%.  In  another  run  the  Ice  surface 
was  cooled  below  the  eutectic  temperature  (-21°C)  causing  solid  salt  to 
crystallize  out  of  the  brine,  forming  a  surface  crust.  The  albedo  was  then 
as  large  as  that  of  dry  snow.  Raising  the  surface  temperature  above  -21®C 
caused  the  salt  to  redissolve,  and  the  albedo  dropped  by  about  30%. 

Growth  rate  dependence  was  tested  by  comparing  Ice  samples  of  the  same 
thickness  and  temperature,  but  which  were  grown  at  different  air  temperatures 
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between  -10®C  and  -37“C.  Spectral  albedos  at  all  wavelengths  were  larger 
for  ice  which  had  grown  faster.  A  difference  of  about  40%  was  observed 
between  the  extreme  cases  at  all  wavelengths. 

Surprisingly,  spectral  albedos  from  the  intermediate  salinity  experiments 
were  the  same  or  slightly  higher  (up  to  20%)  than  those  from  the  corresponding 
high  salinity  (32  °/oo)  cases.  This  was  caused  by  the  formation  of  a  thin 
brine  skim  on  the  surface  which  was  not  present  on  the  intermediate  salinity 
ice.  The  skim  filled  in  surface  irregularities  and  reduced  backscattering  in 
the  uppermost  ice  layers.  This  effect  was  most  pronounced  for  an  air  temperature 
of  -20°C;  at  an  air  temperature  of  -10°C  the  skim  was  barely  noticeable  and 
albedos  for  the  16  °/oo  and  32  o/qo  cases  were  nearly  the  same. 

Spectral  extinction  coefficients  (xj^)  observed  for  growing  young  ice 
were  from  1.5  to  15  times  greater  than  for  the  thick  summer  ice  studied  in 
our  field  experiments.  To  check  these  results  the  ice  was  heated  with  the 
artificial  sun  and  bulk  extinction  coefficients  were  determined  from  changes 
in  the  vertical  temperature  profiles.  Results  from  the  two  methods  agreed 
to  within  5%,  supporting  the  accuracy  of  the  optical  data.  As  with  the 
albedos,  spectral  extinction  coefficients  were  larger  for  higher  growth  rates 
and  decreased  with  increasing  temperature.  At  500  nm,  ranged  from  about 
2  m"^  for  ice  grown  at  an  air  temperature  of  -10®C  to  12  m“^  at  -37°C.  When 
the  air  temperature  was  raised  to  -2°C  and  the  ice  allowed  to  equilibrate, 

<500  decreased  by  as  much  as  8  m"^.  The  changes  in  extinction  with  temperature 
and  growth  rate  appeared  to  be  smooth  except  when  the  upper  layers  of  the  ice 
crossed  the  eutectic  temperature. 

Since  the  absorption  coefficients  of  brine  and  pure  ice  are  essentially 
independent  of  temperature,  the  optical  properties  of  sea  ice  are  governed 


14 


by  the  effectiveness  of  scattering  In  the  Ice.  Enhanced  scattering  results 
In  higher  albedos;  It  also  Increases  the  extinction  coefficient  because  the 
optical  path  length  Is  greatly  Increased  by  multiple  scattering  with  a 
corresponding  Increase  In  the  probability  that  the  light  will  be  absorbed 
within  the  Ice.  Scattering  In  turn  Is  determined  primarily  by  the  platelet 
structure  and  the  distribution  of  brine  between  the  platelets.  The  amount 
of  scattering  Is  Inversely  related  to  the  temperature.  As  brine  pockets  and 
channels  open  up  with  Increasing  temperature,  the  brine/lce  Interfaces 
coalesce  and  become  more  rounded.  If  the  Ice  temperature  drops  below  the 
eutectic  point,  salt  crystals  are  precipitated  out  causing  a  sudden  Increase 
In  the  scattering.  Growth  rate  Is  Important  because  It  determines  the  platelet 
and  crystal  sizes  and  the  Initial  distribution  of  brine  In  the  Ice.  Albedo  Is 
also  Influenced  by  changes  In  surface  conditions.  Flooding  of  the  surface  due 
to  melting  or  to  concentrated  brine  being  squeezed  out  gives  a  smooth  surface 
and  lowers  the  albedo,  and  melting  followed  by  drainage  leaves  a  crumbly 
surface  with  a  higher  albedo. 

From  these  results  It  appears  that  the  extinction  coefficient  of  natural 
sea  Ice  will  vary  with  depth  and  time  since  the  optical  properties  respond 
to  changes  In  the  temperature  and  salinity  of  the  ice.  Thus,  the  optical 
properties  of  sea  Ice  should  be  In  a  continual  state  of  change.  In  the  Arctic 
this  Involves  a  transient  phase  for  the  first  two  years  until  the  salinity 
profile  stabilizes  followed  by  a  cyclic  variation  with  seasonal  environmental 
changes. 

While  the  results  of  the  present  experiments  Isolate  the  effects  which 
Influence  the  optical  properties,  further  work  Is  necessary  to  give  a 
quantitative  connection  between  a  and  k  and  the  Ice  structure.  The  weak 
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f  dependence  of  the  optical  properties  on  water  salinity  should  be  examined  in 

more  detail.  Finally,  the  effects  of  vapor  bubbles  have  not  been  studied 

I 

here,  but  should  be  included  in  any  general  theory  of  the  optical  properties 
of  sea  ice.  Several  of  these  problems  have  been  studied  in  the  scattering 
experiments  described  below. 

The  complete  results  of  this  work  have  been  released  as  Scientific  Report 
No.  17,  and  a  condensed  version  submitted  to  the  Journal  of  Glaciology  for 
publication.  Me  also  plan  to  write  an  additional  paper  describing  the 
theoretical  part  of  the  work. 

Scanning  Photometer.  Construction  and  testing  of  the  scanning  photometer 
has  been  completed.  The  instrument  scans  over  visible  and  infrared  wavelengths 
*  from  400  to  2500  nm  using  a  circular  variable  interference  filter,  together 

^  with  a  set  of  twelve  fixed  band  interference  filters.  The  spectral  resolving 

power  (wavelength/bandpass)  ranges  from  about  25  at  400  nm  to  46  at  1300  nm. 
Thus,  in  addition  to  extending  the  wavelength  region  to  cover  more  than  98% 
of  the  solar  spectrum,  the  new  system  has  almost  20  times  greater  spectral 
resolution  between  800  and  1000  nm  than  do  our  previous  instruments.  The 
radiation  detectors  include  two  lead  sulfide  photoconductors  for  the  infrared 
beyond  1000  nm  and  a  blue-enhanced,  ultra-low-leakage  photodiode  for  visible 
and  near  infrared  wavelengths.  The  system  will  detect  spectral  irradiances 
of  less  than  0.1  pW/cm^/nm,  giving  2-3%  accuracy  in  the  incident  irradiance 
at  2.5  microns  on  a  heavily  overcast  summer  day  in  the  Arctic. 

4 

The  entire  photometer  is  contained  in  an  aluminum  cylinder  20  cm  in 
^  diameter  and  20  cm  in  height  and  weighs  less  than  35  pounds  including  batteries. 

To  eliminate  A.C.  noise  problems,  only  D.C.  power  supplies  which  run  off 
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rechargeable  lead  acid  batteries  are  used.  C-MOS  electronics  have  been 
Installed  where  possible  for  low  power  drain  and  low  temperature  operation. 

The  photometer  Is  Intended  primarily  for  measuring  Incident  spectral 
Irradlance  and  Infrared  spectral  albedos,  but  It  has  a  fiber  optics  attachment 
If  measurements  within  the  Ice  are  required. 

A  field  experiment  was  carried  out  at  Point  Barrow  from  mid-May  until 
the  end  of  June  1979.  Most  of  the  observations  were  made  on  the  seasonal 
sea  Ice  between  1/2  and  2  miles  offshore  near  the  Naval  Arctic  Research 
Laboratory.  The  goals  of  the  experiment  were:  (1)  to  determine  spectral 
albedos  of  different  surface  types  during  the  onset  of  the  melt  season, 

(11)  to  study  Incident  spectral  Irradlance  versus  cloudiness,  and  (111)  to 
do  some  Initial  testing  of  our  methods  for  the  summer  lead  experiment. 

After  some  Initial  adjustments  the  photometer  performed  extremely  well. 

A  broken  ground  wire  was  the  only  serious  problem  encountered.  At  an 
ambient  temperature  of  -10°C  the  Instrument  would  run  for  3  to  4  hours  on 
one  charge  of  the  battery  pack.  Measurements  of  Incident  spectral  Irradlance 
were  digitized  and  recorded  on  magnetic  tape,  which  was  much  faster  than 
hand  logging.  In  order  to  measure  entire  spectra  quickly  before  cloud 
conditions  could  change  appreciably.  The  scans  were  played  back  through  an 
XY  plotter  at  NARL  In  case  the  digitized  data  were  later  lost.  Initial 
calibrations  of  Instrument  sensitivity  were  carried  out  at  the  lab  at  +20®C 
as  Insurance  against  possible  Instrument  damage  In  shipping.  Since  the 
sensitivity  Is  temperature  dependent  and  a  suitable  controlled  temperature 
environment  was  not  available  at  Point  Barrow,  a  complete  set  of  calibrations 
was  carried  out  this  September  In  our  own  cold  rooms. 


Incident  and  reflected  Irradlances  were  obtained  at  40  different  sites, 
ranging  from  dry  wind-packed  snow  to  well  developed  melt  ponds.  Simultaneous 
Kipp  radiometer  measurements  were  taken  both  to  provide  total  albedos  as  a 
check  for  the  spectral  data,  and  to  monitor  fluctuations  In  the  Incident 
radiation  due  to  changes  In  cloud  conditions. 

The  albedo  data  were  recorded  by  hand  at  52  wavelengths  and  were  reduced 
each  evening.  The  results  are  consistent  with  our  previous  data  from  400  to 
800  nm  and  with  the  Kipp  albedos.  Considerable  Improvement  was  achieved  from 
800  to  1000  nm  with  the  higher  resolution  and  good  spectral  detail  was  observed 
at  longer  wavelengths.  Although  some  differences  exist  as  we  expected,  our 
data  for  snow  appear  to  be  consistent  with  reflectances  measured  by  O'Brien 
and  Munis  In  the  laboratory. 

Computer  processing  of  the  results  Is  now  nearly  complete,  and  the  data 
have  been  recovered  from  all  sites  but  one.  From  the  XY  plots  It  Is  evident 
that  clouds  do  absorb  much  of  the  Infrared  radiation;  however,  a  good  deal 
more  infrared  energy  reaches  the  ground  than  we  originally  estimated. 

Reduction  and  Interpretation  of  the  data  will  be  completed  next  year  and 
the  results  will  then  be  written  up  for  publication. 

A  short  talk  on  the  design  and  operation  of  the  photometer  was  presented 
before  the  Northwest  Glaclological  Society  this  November.  This  Information 
will  also  be  written  up  In  a  short  paper  for  the  Instruments  and  Methods 
section  of  the  Journal  of  Glaciology. 

Scattering  Experiments.  Measurements  of  the  light  scattering  properties 
of  sea  ice  have  been  completed.  Fifty-eight  cases  have  been  studied 
Including  bubble-free  Ice,  fresh  bubbly  Ice  (glacier  Ice),  salt  ice  grown  at 
various  temperatures  In  the  cylindrical  tank,  and  grease  Ice. 
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Light  scattering  is  characterized  by  a  coefficient,  o,  defined  as  the 
fraction  of  light  scattered  per  unit  length  out  of  a  collimated  beam  incident 
on  the  test  sample,  and  by  the  phase  function  which  describes  the  angular 
distribution  of  the  scattered  radiation.  These  two  parameters  define  the 
role  of  scattering  in  the  transfer  of  radiation  through  sea  ice. 

For  bubble-free  ice  the  apparatus  originally  detected  faint  scattering 
from  surface  irregularities.  When  these  were  removed  by  polishing  the  sample 
faces,  no  scattered  light  could  be  detected;  thus  the  scattering  coefficient 
was  less  than  10~^  m~^  and  a  much  larger  sample  would  be  needed  to  measure 
the  scattering.  However,  since  the  absorption  coefficient  is  several  orders 
of  magnitude  larger  (more  than  4  x  lO"^  m"^)  than  a,  scattering  can  be  ignored. 

For  sea  ice  the  scattering  is  much  more  efficient;  measured  values  of  a 
lie  between  0.03  and  0.4  m~-.  The  phase  functions  show  strong  forward 
scattering  with  a  minimum  near  120®  and  a  slight  increase  for  larger  deflection 
angles.  The  phase  function  is  independent  of  wavelength  from  400  nm  to  at 
least  800  nm  which  suggests  that  scattering  is  controlled  by  the  refractive 
indices  of  salt  water,  ice,  and  air,  none  of  which  are  strongly  wavelength 
dependent. 

Substantial  variations  in  the  scattering  coefficients  and  phase  functions 
were  observed  for  different  salt  ice  samples.  A  strong  temperature  dependence 
between  -10®C  and  -30®C  was  found  where  scattering  coefficients  for  the  colder 
samples  were  2  to  3  times  larger.  The  phase  functions  showed  some  increase 
in  forward  and  backscattering  with  decreasing  temperature  as  well  as  a  weak 
dependence  on  sample  orientation.  Large  increases  in  o  occurred  when  the 
ice  was  cooled  below  the  eutectic  point  and  the  scattering  became  nearly 
isotropic. 


Scattering  coefficients  for  grease  ice  which  had  been  frozen  solid  were 
somewhat  larger  than  for  samples  from  the  thin  ice  experiments  at  the  same 
temperature  and  showed  more  sidescattering  and  backscattering.  This  is 
reasonable  since  the  grease  ice  is  made  up  of  many  randomly  oriented  small 
crystals  rather  than  a  smaller  number  of  large  columnar  crystals  characteristic 
of  the  interior  of  the  thin  ice. 

Results  thus  indicate  that  the  scattering  coefficient  depends  on  both 
the  freezing  rate  and  the  temperature  of  the  samples.  The  phase  function, 
however,  is  not  strongly  influenced  by  these  quantities.  This  should  simplify 
the  parameterization  of  the  scattering  properties.  Final  analysis  is  in 
progress  and  the  results  are  being  prepared  for  publication  in  the  Journal 


Optical  Properties  of  Cold  Seasonal  Ice  and  Grease  Ice.  The  data  from 
the  spring  field  experiment  carried  out  near  Prudhoe  Bay  have  been  reduced. 
This  study  gave  us  the  opportunity  to  measure  the  optical  properties  of  cold, 
thick  first-year  ice  which  had  not  undergone  a  melt  season.  Spectral  albedos 
between  400  and  1000  nm  for  cold  ice  with  the  snow  scraped  from  the  surface 
were  quite  high,  even  though  no  surface  granular  layer  was  present.  They  lie 
between  our  previous  results  for  melting  first-year  white  ice  and  melting 
multiyear  white  ice.  The  total  albedo  was  0.55  as  compared  with  0.47  and  0.57 
for  first-year  and  multiyear  ice,  respectively. 

Extinction  coefficients  (xj^)  for  cold  ice  were  larger  than  for  the 
interior  of  summer  ice,  but  smaller  than  for  the  young  ice  grown  in  the 
laboratory.  Values  of  at  500  nm  ranged  from  1.3  to  2.0  m“^  as  compared 
with  0.6  to  1.2  m"^  for  summer  ice  and  2.5  to  12  m"^  for  young  ice.  In  one 
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case,  however,  <500  was  comparable  (0.85  m~M  to  the  summer  values.  At  this 
site  the  uppermost  16  cm  consisted  of  nearly  fresh  ice  from  the  nearby 
Shaviovik  River. 

The  present  results  are  consistent  with  the  findings  of  the  thin  ice 
experiments  which  show  that  a  and  k  are  larger  for  colder  more  rapidly  grown 
ice.  Thus,  since  the  average  growth  rate  of  thick  cold  first-year  ice  is 
slower  than  for  thin  young  ice,  its  extinction  coefficients  and  hence  light 
scattering  within  the  ice  should  be  less  than  for  the  young  ice  but  greater 
than  for  melting  summer  ice.  The  corresponding  comparison  for  albedos  is 
not  as  simple  because  variations  in  the  surface  conditions  such  as  flooding 
on  the  formation  of  a  drained  granular  layer  cause  large  changes  in  the 
albedos.  While  a  quantitative  analysis  will  be  required  to  compare  albedos 
for  ice  of  different  thicknesses,  the  present  albedos  also  appear  to  be 
consistent  with  our  previous  results. 

The  laboratory  studies  of  grease  ice  indicate  that  the  albedos  are 
considerably  lower  than  for  the  same  thickness  of  solid  ice  grown  in  the 
cylindrical  tank  at  similar  air  temperatures.  The  extinction  coefficients 
on  the  other  hand  are  only  slightly  smaller.  This  suggests  that,  in  addition 
to  the  influence  of  the  very  high  "brine  volume"  of  grease  ice,  its  light 
scattering  properties  are  significantly  different  from  those  of  ice  grown  at 
the  bottom  of  an  existing  ice  layer.  A  quantitative  investigation  of  the 
comparison  of  cold  young  ice,  grease  ice,  and  thick  summer  ice  will  be  an 
important  part  of  the  parameter  study  planned  for  next  year.  We  intend  to 
use  the  scattering  data  in  conjunction  with  the  16-stream  radiative  transfer 
model  to  determine  how  varying  the  concentrations  of  the  different  scattering 
inhomogeneities  affects  the  albedo  and  transmission. 
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FLUID  MECHANICS  AND  SEA  ICE 

During  the  past  year  graduate  student  Edward  Josberger  completed  his 
thesis  on  the  ablation  of  vertical  ice  walls  submerged  in  warm  seawater. 

S.  Martin  and  P.  Kauffman  completed  an  analysis  of  the  laboratory  studies 
of  the  interaction  of  waves  and  grease  ice,  and  carried  out  a  cooperative 
field  program  in  March  1979  on  the  properties  of  the  Bering  Sea  ice  edge 
with  NOAA  and  the  Scott  Polar  Research  Institute.  Another  graduate  student 
supported  on  the  contract,  Ms.  Jane  Bauer,  began  a  series  of  laboratory 
experiments  on  the  melting  of  ice  floes  in  a  wave  field. 

Our  reports  and  papers  during  the  past  year  include  the  technical  report 
by  Edward  Josberger  entitled  "Laminar  and  Turbulent  Boundary  Layers  Adjacent 
to  Melting  Vertical  Ice  Walls  in  Salt  Water".  This  has  been  revised  by 
Josberger  and  Martin  into  the  paper  entitled  "A  Laboratory  and  Theoretical 
Study  of  the  Boundary  Layer  Adjacent  to  a  Vertical  Melting  Ice  Wall  in  Salt 
Water"  and  submitted  to  the  Journal  of  Fluid  Mechanics.  S.  Martin  and 
P.  Kauffman  have  in  the  final  stages  of  preparation  a  manuscript  entitled 
"A  Field  and  Laboratory  Study  of  Wave  Damping  by  Grease  Ice"  which  will 
shortly  be  submitted  to  the  Journal  of  Glaciolodv.  This  is  an  account  of 
our  laboratory  experiments  and  field  observations  on  wave  damping  by  sea  ice, 
and  gives  as  complete  an  account  as  we  are  able  to  provide  of  the  role  of 
grease  ice  in  the  marginal  seas. 

Our  major  effort  during  the  past  year  was  the  organization  of  an  ice 
edge  cruise  on  the  NOAA  ship  SURVEYOR.  Although  this  work  was  primarily 
supported  by  the  OCSEAP  program,  the  present  ONR  contract  contributed  to 
the  cruise  in  two  ways.  First,  the  contract  paid  for  the  costs  related  to 
the  borrowing,  insuring,  and  repair  of  a  wave-rider  buoy,  which  was  deployed 
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from  the  ship  to  gather  wave  data  in  support  of  a  wave  attenuation  study 
carried  out  by  Dr.  Vernon  Squire  and  Mr.  Stuart  Moore  of  the  Scott  Polar 
Research  Institute.  Squire  and  Moore  have  written  up  part  of  this  Bering 
Sea  work  in  a  paper  entitled  "Direct  Measurement  of  the  Attenuation  of 
Ocean  Waves  by  Pack  Ice",  which  is  in  press  at  Nature. 

Second,  Ms.  Bauer  carried  out  a  study  on  the  cruise  of  the  relative 
motion  of  the  bands  of  ice  which  form  at  the  ice  edge.  These  bands  measure 
about  1  km  wide  by  10  km  long,  and  form  during  periods  of  off-ice  winds. 

In  Bauer's  study  she  mounted  six  visual  targets  on  different  ice  floes 
within  1  km  of  the  edge  following  a  period  of  easterly  winds  which  compacted 
the  ice  edge.  When  the  winds  shifted  to  the  northeast,  the  ice  moved  to  the 
southwest  with  ice  bands  starting  to  form.  She  then  tracked  the  motion  of 
the  targets  from  helicopter  overflights  with  the  Global  Navigation  System 
over  a  23  hour  period.  Her  results  showed  that  the  band  moved  to  the 
southwest  at  about  of  the  wind  speed  and  at  35®  to  the  right  of  the  wind. 

At  the  beginning  of  her  observations  the  ice  in  the  bands  was  part  of  the 
pack;  after  23  hours  the  ice  was  at  least  25  km  to  the  southwest  of  the 
pack. 

Our  observations  show  that  these  bands  form  because  the  ice  floes  near 
the  edge  have  a  greater  aerodynamic  roughness  than  the  interior  ice.  The 
cause  of  the  roughness  increase  is  that  the  floes  near  the  edge  are  rafted 
and  ridged  by  the  action  of  ocean  swell.  When  the  winds  blow  off  the  ice,  the 
rough  ice  moves  south  faster  than  the  interior  ice,  which  leads  to  the 
formation  of  the  bands.  Bauer  and  Martin  include  a  discussion  of  this  work 
in  a  chapter  entitled  "Bering  Sea  Ice  Edge  Phenomena",  which  is  written  for 
a  proposed  Bureau  of  Land  Management  book.  We  plan  to  rewrite  this  chapter 
in  a  more  technical  manner  for  submission  to  the  Journal  of  Geophysical  Research. 
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Finally,  we  completed  construction  of  our  new  1  m  deep  by  1  m  wide 
by  5  m  long  wave  tank.  Using  this  tank,  Ms.  Bauer  has  begun  a  series  of 
experiments  on  how  Ice  floes  melt  when  they  are  placed  In  warm  salt  water 
through  which  waves  propagate.  The  purpose  of  this  experiment  Is  to 
discover  how  the  wave- Induced  boundary  layer  under  the  Ice  floe  enhances 
the  Ice  ablation,  then  to  compare  the  experimental  results  with  theory. 
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1.  Grenfell,  T.  C.,  The  effects  of  ice  thickness  on  the  exchange  of  solar 
radiation  over  the  polar  oceans.  Journal  of  Glaciology,  22,  87, 
305-20,  1979. 


Total  transmission,  absorption,  and  reflection  of  solar  radiation  have 
been  determined  for  bare  blue  and  white  ice  between  0.02  and  0.8  m  in 
thickness  as  well  as  for  blue  ice  covered  with  0.01  to  0.4  m  of  dry 
packed  snow.  The  calculations  were  performed  at  45  wavelengths  between 
400  nm  and  2150  nm  using  a  two-stream  model  to  account  for  the  finite 
thickness  of  the  ice  and  snow  layers.  Total  radiative  energies  were 
found  by  numerical  integration  over  wavelength.  The  results  were 
compared  with  corresponding  calculations  for  optically  thick  ice  of  the 
same  types.  Albedos  increase  from  about  0.05  for  open  water  to  a 
maximum  of  0.9  for  thick  snow.  For  0.8  m  blue  and  white  ice,  predicted 
albedos  on  cloudy  d6(ys  are  0.28  and  0.67  respectively.  Under  clear 
skies  these  albedos  decrease  by  10  to  30%.  Total  transmission  through 
thin  ice  (less  than  0.8  m)  is  from  50%  to  300%  greater  than  is  predicted 
by  Beer's  law  depending  on  ice  type  and  cloud  cover.  Radiative  energy 
absorption  at  the  surface  is  independent  of  thickness,  but  significant 
departures  from  Beer's  law  of  as  much  as  200%  are  evident  in  all  cases 
below  a  depth  of  2.5  mm.  A  parameterization  scheme  is  presented  for 
incorporating  these  results  into  heat-  and  mass-balance  studies. 


2.  Niedrauer,  T.  M.,  and  S.  Martin,  An  experimental  study  of  brine  drainage 
and  convection  in  young  sea  ice.  Journal  of  Geophysical  Research,  84, 
C3,  1176-86,  1979. 


In  a  series  of  experiments  using  a  1 .6-tnn-thick  freezing  tank,  thin 
sections  of  salt  water  ice  were  grown  which  exhibit  the  same  drainage 
features  as  natural  sea  ice.  The  tank  design  permitted  photographs  to 
be  taken,  while  thermocouples  mounted  in  the  tank  walls  recorded  the 
temperature  profiles  within  the  ice.  Convection  was  observed  in  both 
*  the  skeletal  layer  and  in  the  brine  channels  by  the  flow  of  dyed  brine. 

Flow  in  the  skeletal  layer  was  cusplike  in  appearance,  consisting  of 
narrow  downflow  regions  separated  by  broad  upflow  regions.  Above  the 
(  skeletal  layer,  several  brine  channels  were  also  usually  present  in  the 

ice,  and  convective  overturning  occurred  in  these  channels.  The 


convection  caused  temperature  fluctuations  of  0.05®C,  which  calculations 
show  Increase  the  vertical  heat  flux  by  1%.  The  brine  drainage  channels, 
which  were  usually  sloped  30°  to  60°  to  the  horizontal,  always  had 
Isotherms  tilted  from  0°  to  13°  In  the  same  direction.  The  brine  channels 
move  both  horizontally  and  vertically  through  the  Ice  by  melting  their 
lower  walls  and  freezing  on  the  upper  walls.  An  analysis  based  on  the 
heat  flux  due  to  brine  channel  convection  shows  that  convection  can  drive 
these  wall  movements.  Our  observations  suggest  that  most  of  the  brine 
movement  In  the  channels  Is  caused  by  recirculation  of  water  from  below 
the  ice.  On  a  small  scale  we  also  observed  the  formation  of  brine  pockets 
from  brine  tubes. 


Josberger,  E.  G.,  Laminar  and  turbulent  boundary  layers  adjacent  to  melting 
vertical  Ice  walls  In  salt  water.  Scientific  Report  No.  16,  Office  of 
Naval  Research,  Contract  N00014-76-C-0234,  Department  of  Atmospheric 
Sciences,  University  of  Washington,  Seattle,  Washington,  185  pp.,  1979. 


A  laboratory  study  of  melting  vertical  ice  sheets  Itnnersed  In  sodium 
chloride  solutions  of  oceanic  temperatures  and  salinities  found  that  the 
convective  motions  generated  by  the  melting  consist  of  two  regions 
confined  to  a  boundary  layer  adjacent  to  the  ice.  The  first  region,  at 
the  bottom  of  the  Ice,  Is  laminar  and  bidirectional,  with  a  2  to  3  mn 
thick  upward  flowing  layer  next  to  the  Ice  Inside  of  a  10  nrni  thick  outer 
downward  flowing  layer.  This  bidirectionality  results  from  the  large 
difference  between  the  thermal  and  saline  diffusivltles  which  contain 
dilute  water  near  the  Ice  and  allows  the  cooling  to  diffuse  further  from 
the  Ice.  Near  the  Ice,  dilute  cold  water  rises;  while  away  from  the  ice, 
cold  saline  water  sinks.  Second,  further  up  the  Ice,  the  Inner  layer 
goes  through  a  region  of  transition  and  becomes  fully  turbulent  when  the 
saline  Grashof  number  reaches  2  x  10®.  The  turbulent  flow  Is  upward, 
unidirectional  and  the  dominant  flow  In  the  oceans. 

Ice-water  Interface  temperature  measurements  over  the  Ice  surface  show  a 
vertically  varying  temperature  In  the  laminar  region  and  a  uniform 
temperature  In  the  turbulent  region.  The  measured  melt  rates  were  highest 
In  the  transition  region,  lowest  In  the  laminar  region,  and  had  a  slow 
vertical  variation  In  the  turbulent  region.  For  water  temperatures 
greater  than  25°C,  the  turbulent  flow  reversed  direction  and  a  laminar 
bidirectional  region  formed  at  the  top  of  the  Ice. 

In  a  theoretical  study  of  the  turbulent  flow,  an  eddy  diffusivlty  models 
the  turbulent  transport  of  momentum,  salt  and  heat  In  the  conservation 
equations  and  the  laboratory  measurements  determine  the  model  parameters. 
The  turbulent  model  consists  of  three  regions:  a  laminar  Inner  region. 
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an  Intermediate  region  where  the  eddy  dlffusivity  increases  linearly 
with  distance  from  the  ice  and  an  outer  region  of  constant  eddy 
dlffusivity.  A  similarity  transformation  reduces  the  governing  partial 
differential  equations  to  a  system  of  ordinary  differential  equations 
which  are  then  numerically  integrated. 

The  numerical  results  show  the  following:  The  eddy  dlffusivity  depends 
parabolically  on  the  density  difference  across  the  boundary  layer.  The 
major  effect  of  the  turbulence  is  to  Increase  the  transport  of  dilute 
water  away  from  the  interface.  The  dilution  effect  overwhelms  the 
thermal  effect  on  density  to  produce  a  unidirectional  upward  flow. 
Finally,  the  injection  of  melt  water  into  the  boundary  layer  reduces  the 
salt  flux  to  the  ice  which  lowers  the  interface  salinity  as  the  water 
temperature  increases.  Use  of  the  results  to  form  a  zero  net  buoyancy 
condition  on  the  buoyancy  integrated  across  the  boundary  layer  predicts 
the  reversal  of  the  turbulent  flow  which  accurately  fits  all  of  the  known 
data. 

Salinity  and  temperature  measurements  with  depth  in  the  vicinity  of  an 
iceberg  in  the  North  Atlantic  Ocean  show  variations  in  these  quantities 
near  the  iceberg  that  may  be  the  result  of  the  upward  flowing  turbulent 
boundary  layer. 


4.  Perovich,  D.  K.,  The  optical  properties  of  young  sea  ice.  Scientific 
Report  No.  17,  Office  of  Naval  Research,  Contract  N00014-76-C-0234, 
Department  of  Atmospheric  Sciences,  University  of  Washington,  Seattle, 
Washington,  151  pp.,  1979. 


Eight  laboratory  experiments  were  performed  to  determine  the  optical 
properties  of  young  salt  ice  and  to  examine  correlations  between  the 
optical  properties  and  the  state  of  the  ice.  Ice  grown  at  different 
temperatures  (-10,  -20,  -30,  and  -37“C)  and  from  water  of  different 
salinities  (0,  16,  and  31  °/oo)  was  investigated.  The  experiments  were 
conducted  using  a  tank  system  which  was  designed  to  grow  ice  similar  to 
that  found  in  nature,  allow  an  accurate  determi nation  of  the  state  and 
structure  of  the  ice,  and  permit  in  situ  optical  measurements  to  be  made. 
Measurements  of  incident,  reflecteJ,  and  transmitted  irradiances  were 
used  in  conjunction  with  a  modified  Dunkle  and  Bevans  photometric  model 
to  determine  spectral  albedos  and  extinction  coefficients.  The  thin  cold 
ice  of  these  experiments  had  albedos  which  were  comparable  to  the  values 
for  the  thicker  warmer  ice  examined  by  previous  researchers;  however, 
extinction  coefficients  were  1.5  to  15  times  greater.  Qualitative 
relationships  between  the  optical  properties  and  the  physical  state  of 
the  ice  were  observed.  As  the  ice  temperature  decreased  (and  the  brine 
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volume  decreased)  both  albedo  and  extinction  coefficient  Increased; 
when  the  ice  temperature  dropped  below  the  eutectic  point  the  Increase 
was  drastic.  The  dependence  of  albedo  and  extinction  coefficient  on 
the  brine  content  of  the  ice  was  found  to  be  complex  with  both  the  brine 
volume  and  its  vertical  distribution  being  significant.  Variations  in 
ice  salinity  over  the  range  4  o/qq  to  14  o/qo  did  not  influence  the 
optical  properties.  Increases  in  albedo  and  extinction  coefficient  were 
primarily  a  result  of  changes  in  the  ice  structure  which  enhanced 
scattering. 

A  four  stream  discrete  ordinates  photometric  model  was  developed  to 
treat  the  case  of  a  floating  ice  slab.  The  model  included  both  aniso¬ 
tropic  scattering  and  refraction  at  the  boundaries.  The  effects  on 
albedo  and  transmittance  of  variations  in  such  model  parameters  as  the 
single  scattering  albedo  and  the  phase  function  were  investigated. 

Using  one  and  two  layer  models,  theoretical  albedos  and  transmittances 
were  compared  to  experimental  values.  The  four  stream  model  was  further 
applied  to  investigate  such  topics  as  the  depth  dependence  of  irradiance, 
cases  of  direct  beam  incident  radiation,  the  spectral  dependence  of 
albedo  and  transmittance,  and  the  effects  of  surface  layers.  It  was 
suggested  that  future  laboratory  research  be  oriented  towards  determina¬ 
tion  of  single  scattering  albedos  and  phase  functions  for  an  extensive 
range  of  ice  types. 
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LIBRARY  CODE  2124 
MONTEREY.  CA  93940 

OR  ARTHUR  LACHENBRUCH 
BRANCH  OF  GEOPHYSICS 
US  GEOLOGICAL  SURVEY 
345  MIDDLEFIELD  ROAD 
MENLO  PARK.  CA  94025 

DR  W  M  SACKINGER 
DEPT  OF  ELECTRICAL  ENGINEERING 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

DR  KEITH  MATHER 
705  GRUENING  BUILDING 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

OR  CHARLES  E  BEHLKE,  DIRECTOR 
INST  OF  ARCTIC  ENVIRONMENTAL  ENG 
UNIVERSITY  OF  ALASKA 
COLLEGE,  AK  99701 

DR  DONALD  W  HOOD 
INST  FOR  MARINE  SCIENCE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
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DR  GUNTER  WELLER 
GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
COLLEGE,  AK  99701 

LIBRARIAN 

NAVAL  ARCTIC  RESEARCH  LABORATORY 
BARROW,  AK  99723 

MGR  INUVIK  RESEARCH  LABORATORY 
BOX  1430 

INUVIK,  NORTHWEST  TERRITORIES 
XOE  OTO  CANADA 

DR  SVENN  ORVIG 
MCGILL  UNIVERSITY 
DEPT  OF  METEOROLOGY 
PO  BOX  6070 

MONTREAL  101,  PQ,  CANADA 

MARINE  SCIENCE  CENTRE  LIBRARY 
MCGILL  UNIVERSITY 
PO  BOX  6070 

MONTREAL  101,  PQ,  CANADA 

DR  E  R  POUNDER 
RUTHERFORD  PHYSICS  BUILDING 
MCGILL  UNIVERSITY 
3600  UNIVERSITY  STREET 
MONTREAL,  PQ  H3A  2T8,  CANADA 

MRS  GAIL  NORWOOD 
METEOROLOGY  LIBRARY 
MCGILL  UNIVERSITY 
MONTREAL  101,  PQ,  CANADA 

DR  WESTON  BLAKE  JR 
GEOLOGICAL  SURVEY  OF  CANADA 
DEPT  OF  ENERGY  MINES  RESOURCES 
601  BOOTH  STREET 
OHAWA,  ONTARIO  K1A  0E8  CANADA 

DEPARTMENTAL  LIBRARY-SERIALS 
DEPT  OF  THE  ENVIRONMENT 
OTTAWA,  ONTARIO  K1A  0H3  CANADA 

DEFENCE  RESEARCH  BOARD 
DEPT  OF  NATIONAL  DEFENCE 
190  O'CONNOR  STREET 
OHAWA,  ONTARIO  KIA  0Z3,  CANADA 


MISS  MOIRA  DUNBAR 
DEFENCE  RESEARCH  ESTABLISHMENT 
OHAWA,  NATL  DEFENCE  HEADQUARTERS 
OHAWA,  ONTARIO  KIA  0Z4,  CANADA 

DR  ROBERT  L  RAUSCH 
DEPT  OF  MICROBIOLOGY 
W  COLLEGE  OF  VET  MEDICINE 
UNIVERSITY  OF  SASKATCHEWAN 
SASKATOON,  SASK  S7N  OWO,  CANADA 

DR  E  L  LEWIS 

FROZEN  SEA  RESEARCH  GROUP 
INSTITUTE  OF  OCEAN  SCIENCES 
9860  W  SAANICH  ROAD,  PO  BOX  6000 
SIDNEY,  BC  V8L  462,  CANADA 

DR  A  R  MILNE 

DEFENCE  RES  ESTABLISHMENT  PAC 
FLEET  MAIL  OFFICE 
VICTORIA,  BC,  CANADA 

METEOROLOGICAL  OFFICE  LIBRARY 
LONDON  ROAD 

BRACKNELL,  BERKSHIRE,  ENGLAND 

DR  T  E  ARMSTRONG 

SCOTT  POLAR  RESEARCH  INSTITUTE 

CAMBRIDGE,  CB  2  1ER,  ENGLAND 

THE  LIBRARIAN 

scon  POLAR  RESEARCH  INSTITUTE 
CAMBRIDGE  CB2  1ER,  ENGLAND 

DR  MICHAEL  KELLY 

SCHOOL  OF  ENVIRONMENTAL  SCIENCES 
UNIVERSITY  OF  EAST  ANGLIA 
NORWICH  NR4  7TJ,  ENGLAND 

NATIONAL  INST  OF  OCEANOGRAPHY 
WORMLEY,  GOOALMING 
SURREY,  ENGLAND 

CENTRE  NATL  DE  LA  RECHERCHE  SCI 
LABORATOIRE  DE  GLACIOLOGIE 
UNI VERS ITE  DE  GRENOBLE  1 
SERV  DE  GEOPHYS,  2,  RUE  TRES 
CLOITRES  38-GRENOBLE,  FRANCE 
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CENTRE  O' ETUDES  ARCTIQUES 
ECOLE  PARTIQUE  DES  HAUTES  ETUDES 
VI  SECTION  (SORBONNE) 

6.  RUE  DE  TOURNON.  PARIS  -6 
FRANCE 

DR  KOU  KUSUNOKI 
POLAR  RESEARCH  CENTER 
NATIONAL  SCIENCE  MUSEUM 
KAGA  1-9-10.  ITABASHI-KU 
TOKYO.  JAPAN 

NORSK  POLARINSTITUTT 
ROLFSTANGVN  12.  POSTBOKS  158 
1330  OSLO  LUFTHAVN 
NORWAY 

PROF  DR  ALEXSANDER  KOSIBA 
KATEFRA  I  OBSERWATORIUM 
METEOROLOGI  I  KLIMATOLOGII 
UNIUERSYTETN  WROCLAWSKIEGO 
WROCLAW  9.  U  CMENTARNA  8 
POLAND 

STOCKHOLMS  UNIVERSITET 
NATURGEOGRAFISKA  INSTITUTIONEN 
BOX  6801 

113  86  STOCKHOLM.  SWEDEN 

DR  ING  0  MAGGIOLO 
FACULTAD  DEINGENIERIA 
HERRERA  Y  REISSIG  585 
MONTEVIDIO,  URUGUAY 


UNIVERSITY  ARCHIVES 
UNIVERSITY  RECORDS  SECTION  BL-10 
UNIVERSITY  OF  WASHINGTON 
SEATTLE.  WA  98195 

RESIDENT  REPRESENTATIVE 
OFFICE  OF  NAVAL  RESEARCH 
ROOM  422.  UNIV  DISTRICT  BLDG  JD-27 
UNIVERSITY  OF  WASHINGTON 
SEAHLE.  WA  98105 

DR  LAWRENCE  COACHMAN 
DEPARTMENT  OF  OCEANOGRAPHY  WB-10 
UNIVERSITY  OF  WASHINGTON 
SEAHLE.  WA  98195 

DR  T  SAUNDERS  ENGLISH 
DEPARTMENT  OF  OCEANOGRAPHY  WB-10 
UNIVERSITY  OF  WASHINGTON 
SEATTLE.  WA  98195 

OR  KNUT  AAGARD 

DEPARTMENT  OF  OCEANOGRAPHY  WB-10 
UNIVERSITY  OF  WASHINGTON 
SEAHLE.  WA  98195 


LIBRARY-CENTRE  FOR  NORTHERN 
STUDIES  AND  RESEARCH 
MCGILL  UNIVERSITY 
1020  PINE  AVENUE  WEST 
MONTREAL.  PQ  H3A  1A2 
CANADA 


